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Abstract 

The phenomenon of particle creation within an almost resonantly vibrating cavity with losses 
is investigated for the example of a massless scalar field at finite temperature. A leaky cavity is 
designed via the insertion of a dispersive mirror into a larger ideal cavity (the reservoir). In the case 
of parametric resonance the rotating wave approximation allows for the construction of an effective 
Hamiltonian. The number of produced particles is then calculated using response theory as well as a 
non-perturbative approach. In addition we study the associated master equation and briefly discuss 
the effects of detuning. The exponential growth of the particle numbers and the strong enhancement 
at finite temperatures found earlier for ideal cavities turn out to be essentially preserved. The 
relevance of the results for experimental tests of quantum radiation via the dynamical Casimir 
effect is addressed. Furthermore the generalization to the electromagnetic field is outlined. 
PACS: 42.50.Lc, 03.70.+k, ll.10.Ef, ll.10.Wx. 
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I. INTRODUCTION 



Since the pioneering work of Casimir [|lj the phenomena of quantum field theory under 
the influence of external conditions have attracted the interest of many authors, see e.g. [Q. 
The original prediction by Casimir, i.e., the attractive force generated between two perfectly 
conducting objects placed in the vacuum, has been verified in different experimental setups 
with relatively high precision [[|. However, its dynamic counterpart with non- stationary 
boundary conditions inducing interesting effects like the creation of particles out of the vac- 
uum has not yet been observed rigorously in a corresponding experiment. The observation of 
quantum radiation could provide a substantial test of the foundations of quantum field the- 
ory and thus be of special relevance. Generally we understand the term quantum radiation 
to denote the conversion of virtual quantum fluctuations into real particles due to external 
disturbances. For the special case of the external disturbances being moving mirrors this 
phenomenon is known as the Dynamical Casimir effect. 

These striking effects have been investigated by several authors, for an overview see 
e.g. [0, |3|] and references therein. We will focus on the effect of particle creation within a 
constructed - resonantly vibrating - leaky cavity. This case is of special importance for an 
experimental verification of the dynamical Casimir effect since the generation of particles is 
enhanced drastically by resonance effects. Employing different methods and approaches it 
has already been shown for ideal cavities (see e.g. ||) that under resonance conditions (i.e., 
when one of the boundaries performs harmonic oscillations at twice the frequency of one of 
the eigenmodes of the cavity) the phenomenon of parametric resonance (see e.g. ||) will 
occur. In the case of an ideal cavity (i.e., one with perfectly reflecting mirrors) this is known 
to lead to an exponential growth of the resonance mode particle occupation numbers, cf. 

1, 0,n,n,ii. 

In view of this prediction an experimental observation of quantum radiation using the 
dynamical Casimir effect appears to be rather simple - provided the cavity is vibrating at 
resonance for a sufficiently long period of time. However, this point of view is too naive 
since neither ideal cavities do exist nor is it possible to match the external frequency to 
the fundamental eigenfrequency of the cavity with arbitrary precision. Consequently, it is 
essential to include effects of leaks as well as effects of detuning, see also M. 



Investigations concerning effects of losses have been performed for example in [ 10 1 in 1+1 
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space-time dimensions based on conformal mapping methods as developed in ||11|| . However, 
these considerations are a priori restricted to 1+1 dimensions and can not be obviously 
generalized to higher dimensions. In 3+1 dimensions the character of the mechanism gener- 
ating quantum radiation - e.g., the resonance conditions - differs drastically from the 1+1 
dimensional situation. 

More realistic (3+1 dimensional) cavities were considered in \T2\ where the effects of 
losses were taken into account by virtue of a master equation ansatz. However, this master 
equation had not been derived starting from first principles. It has already been noted in 
fill that the employed ansatz is adequate for a stationary cavity - but not necessarily for 
a dynamic one. In addition, most papers did not include temperature effects - which may 
contribute significantly in an experiment. It has been shown in Ref. that for an ideal 
cavity the effect of particle production at finite temperature is enhanced by several orders 
of magnitude in comparison with the pure vacuum contribution. 

In this article we will adopt the canonical approach which has proven to be general, 
successful, and is - in addition - also capable of including temperature effects. However, the 
aforementioned approach still lacks a generalization for systems with losses. We are aiming 



at providing a remedy in this field [ [l"3|j . 

This paper is organized as follows: In Sec. II we present a model system and derive the 
effective Hamiltonian for the resonance case. In Sec. [ITI| we will calculate the number of 
created particles in the cavity after one of the walls has performed resonant oscillations by 



means of response theory. In Sec. [TV] we will derive and solve the associated master equation 
and show consistency with the results obtained in Sec. [Hi]. In Sec. |V] a non-perturbative 
approach is presented and compared with the other results. We derive a treshold condition 



- valid for leaky cavities - for a possible detuning from the fundamental resonance in Sec. |V] 
. We shall close with a summary, a discussion, a conclusion and an outlook. 
Throughout this paper natural units given by h = c = ks = 1 will be used. 
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FIG. 1: Model of a leaky cavity. A large ideal cavity is split up by a dispersive mirror into a leaky 
cavity and a reservoir. The left (ideal) wall of the cavity is vibrating. 

II. GENERAL FORMALISM 
A. The Leaky Cavity 

We want to investigate the effects of a non-ideal cavity in view of the dynamical Casimir 
effect. For that purpose we have to construct a suitable model system. One simple way 
to do that is to insert a dispersive mirror into an ideal cavity while keeping all other walls 
perfectly reflecting. Thereby two leaky cavities are formed. Particles in the left imperfect 
cavity are now able to leave into the right larger box (the reservoir) . For reasons of simplicity 
we consider a rectangular cavity as depicted in Fig. |]. The setup in Fig. [I] is not a new 
idea. A similar - but static - system has already been treated in [111, fll|. However, here in 



addition the left wall is moving with a prescribed trajectory during the time interval [0, T]. 
For ideal cavities this is known to lead to a squeezing of the vacuum state which causes the 
creation of particles inside the cavity, see e.g. 0. 

Note that we are assuming a finite reservoir with a discrete spectrum instead of an infinite 
one leading to a continuum of modes. Since, in an experimental setup, the vibrating cavity 
will most likely be surrounded by walls, etc. (imposing additional boundary conditions), 
this assumption should be justified. 

Assuming a surrounding perfectly reflecting wall is a first idealization of the real situation. 
However, in order to minimize the error obtained by this procedure the experiment could 
be designed in this way, see also Fig. § in section [IX| below. 

The ideal mirrors can be simulated by infinitely high potential walls inducing Dirichlet 
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V(x) 
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FIG. 2: Diagram of the x-dependence of the used potential. 

boundary conditions. For the additional dispersive mirror we use the 5-type model potential 
proposed in Ref. 0, [17 ] 



V(x;t) 



jS(x — b) if a{t) < x < c 
oo otherwise 



(1) 



see also Fig. ^. The parameter 7 represents the transmittance of the internal mirror, whose 



reflection and transmission amplitudes are determined as [|16| 



?7 



r 



00 



(2) 



io + 27 io + 27 

Note that the general procedure presented in this article is independent of the particular 
form of the potential - the aforementioned one has just been chosen for convenience. For 
a more realistic scenario one could apply square-well or Gaussian potentials. In a realistic 
experiment where one would want to create photons instead of scalar particles a dispersive 
mirror could be realized using a thin dielectric slab with a very high dielectric constant. Such 
a mirror could then be approximated by a space-dependent permittivity e(x) = 1 + aS(x). 
This will be addressed in Sec. |XB| . 



B. Hamiltonian 



Throughout this article we will use the notation of [|18j where the particle production in 
an ideal vibrating cavity was calculated - for a more general treatment see e.g. |19[]. We 
consider a massless and neutral scalar field coupled to an external potential: 

1 



C = -(d^)(d^)-V^ 2 . 



(3) 



The perfect mirrors can be incorporated by imposing the corresponding boundary condition 
on $. By expanding the field $ 

$(r,t) = Y,Q»(t)Mr;t) (4) 

into a complete and orthonormal set of functions / M (r; t) satisfying 

J d 3 rf;(r,t)f u (r,t) = S, u , (5) 
J2f;(r;t)U(r';t) = 5 3 (r-r>), (6) 

{2V-A}/ M (r;t) = ^(t)/ |l (r;t) (7) 

one can reach a more convenient form suitable for doing calculations. Since $ is a real 
field, we can choose the set / M to be real. Note that the time dependence of eigenfunctions 
and eigenfrequency is solely induced by the moving boundary. Inserting this expansion into 
Eq. @ transforms the Lagrangian into ]I5| 



H fj, 

fJ,U 

+- Q»M^{t)M VR {t)Q v , (8) 



2 

/J.UK 



where M^ u (t) is an antisymmetric matrix given by 



M, u = J dh^{T-t)f v {r-t). (9) 

This matrix describes the coupling strength between two different modes. We introduce the 
canonical conjugated momenta 

BL 

Furthermore we apply the usual Legendre transform to a Hamiltonian representation and 
perform the quantization. This yields 

k = \ E p l + \ E + E P » M ^» ■ ( n ) 



6 



The above Hamiltonian can be sub-classified into 



H = H + Hj + Hj , 



(12) 



where the single Hamiltonians are given by 



(13) 




±5>n£(f)Q£, 



(14) 




(15) 



The deviation AQ^(t) = D^(t) — (fi°) 2 denotes the difference of the (squared) time-dependent 
eigenfrequencies f^(i) from the unperturbed ones (fi°) 2 . The first term H is the Hamil- 
tonian of harmonic oscillators. The remaining terms will further on be called squeezing 
interaction Hamiltonian and velocity interaction Hamiltonian. We want to point out that in 
the case of a static system (where the eigenfunctions and eigenfrequencies fl^ are constant 
in time) the complete interaction Hamiltonian Hi = Hf + Hj will vanish. The derivation 
of the eigenfunctions f^(r,t) and M^ v {t) will be treated in the following subsection. 

C. Eigenmodes 

As has already been mentioned, we want to find a set of functions satisfying {2V(r; t) — 
A}/ M (r; t) = r2 2 (t)/ M (r; t). Any time dependence can only be induced by the moving bound- 
aries. At first we will just consider the spatial dependence i.e., the stationary problem. The 
differential equation can be treated using the separation ansatz / M (r ) = f^(r^)f^(r±) where 
fn(r\\) depends only on the coordinate parallel to the wall velocity and f tt (r±) is dependent 
on the perpendicular coordinates. For the special case of our model system this means 
fn( r ) — fu( x )fji(y)fu( z ) leading to the trivial y and z dependence of the eigenfunctions 



with Ay and Az denoting the dimensions of the cavity and the frequencies relating via 




(16) 



(17) 




(18) 
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The remaining differential equation reads 

{2 7 5(a; - b) - dl) /» = (^) 2 /» (19) 

where the Dirichlet boundary conditions coming from the perfect mirrors on either side can 
be satisfied by the ansatz 

L M sin[f2^(x — a)] if a < x < b 
fjl(x) = \ R P sin[0*(c - x)} if b < x < c ■ (20) 
elsewhere 



The eigenfunctions have to obey the continuity conditions 



ffa I b) - f*(x 1 b) = 0, (21) 

QfX QfX 

^i fe )-^T&) = 2 7 #(&), (22) 



where the latter can be obtained via integration. These conditions can be combined to an 
eigenvalue equation for fl^ 

- %L = cot [Oftb - a)} + cot [0*(c - b)} = -1 . (23) 

Though there is no obvious analytical solution of this equation, a numerical solution can 
always be obtained for given cavity parameters {0,6,0,7}. However, via introducing the 
dimensionless perturbation parameter 77 = it is also possible to obtain an approximate 
analytical solution. Note that this parameter is small 77^ 1 in the limit of the internal 
mirror being nearly perfectly reflecting. Since the trigonometric functions are very sensitive 
to small frequency variations one can solve the equation using a series expansion in r? M . It 
is obvious that if the right hand side goes to —00 one of the addends or even both can 
become relevant. This depends on the ratio (b — a)/(c — b) and its inverse which are both 
assumed to be non-integer numbers in the following non-perturbative calculations implying 
that only one of the addends is dominating. Accordingly, expanding around the poles of one 
addend one yields a polynomial that can be solved for as a series expansion in rj^ 1. 
Depending on the chosen addend one obtains two sets of approximate eigenfrequencies 

n , _ JVL _ 1 n 



b — a 2(6 — a) 

1 / c — b\ „ 

— r cot n x Tt- r] n . 

4(6 -a) V b-aj ln *> 1 



Q x = J^. !_„ 

n - r c-b 2(c-6) /r ' 
1 / 6 



co t I n x 7T— — j ry* 



4(c-6) V * c 

+ ^ «,r) » (24) 

which constitute a determining polynomial for Note that the index \i = (n x ,l/r) is 
a multi-index, where / and r stand for left-dominated and right-dominated, respectively. 
However, it can be shown easily that the quality of the linear (in rj) approximation suffices 
already for moderate values of 7 > 50. The insertion of (p4|) into the ansatz (p0|) leads to two 
classes of eigenfunctions: left-dominated and right-dominated, respectively. The differences 
between those are clearly visible in Fig. 0. 

In order to avoid the confusion arising from a set of perturbation parameters {rj^} we 
will introduce the fundamental one via 

1? = *, = ^, (25) 
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to which all others are evidently related via rj^ = Q^/Q^rj. Note that this distinction between 
the classes of eigenfunctions is applicable only for small values of r\. 

Consequently, the eigenfunctions can be labeled by multi-indices fi = (n x ,n y ,n z ,r/l): 3 
quantum numbers n xyz G N + and a flag r/l denoting the class (right- or left-dominated, 
respectively) of the eigenfunction. Now we want to consider the effect of one moving bound- 
ary. It is taken into account by substituting a — > a(t) everywhere in the eigenmodes and 
-frequencies. Thereby a time dependence of the eigenfunctions as well as of the eigenfre- 
quencies is introduced. This induces a non- vanishing coupling matrix M^(t) as well as the 
frequency deviation AQ^(t). For small oscillations of the boundary 

a(t) = a + e(b — a ) sin(a>£) (26) 

with a small amplitude e 1 it will be useful to separate the time dependence using 

M^(t) = a(t) J d 3 r^(r;t)Ur;t) 

= a(t)m^(t). (27) 

The geometry factor m^ v {t) is approximately constant m^ v {t) = m^ u + O (e) in this case. 
Consequently, one is lead to 

M^{t) = m^a(t) + O (e 2 ) . (28) 
9 




FIG. 3: Illustration of the lowest left- and right-dominated eigenmodes ff r n{x) for 771^ = 0.1. 



Since the time-dependence of the right-dominated modes is less complicated than that of 
the left-dominated ones, it is advantageous to exploit the antisymmetry of M^ v which also 
implies an antisymmetry of m^ v . For the following calculations the coupling of the lowest 
left-dominated mode \i = (1,1,1,/) to some right-dominated one v = (n x ,n y ,n z ,r) will 
be of special relevance. The y and z integrations simply generate Kronecker symbols and 
therefore the geometry factor results as 

df x 



" da 



-8l, n Jl,n, / dxf l 
J a 



c-b n*, 



b) sin I n x n 



c-b 



+0 (v 2 ) 



(29) 



D. Canonical Quantization 



Aiming at the calculation of possible particle creation effects (expectation values of par- 
ticle number operators) it is convenient to introduce the creation and annihilation operators 

1 



fijQ M (t)+*P M (t) 



(30) 
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obeying the usual bosonic equal time commutation relations 

[al(t),al(t)] = 

These operators diagonalize the free Hamiltonian 



0. 
0. 



(31) 



(32) 

The following calculations will most conveniently be done in the interaction picture where 
the dynamics of an observable is governed by Hq 



dY 
~dt 



H ,Y 



+ 



dY 
~dt 



(33) 



explicit 



For reasons of generality and to include finite temperature effects we describe the state 
of a quantum system by a statistical operator whose dynamics is determined by the von 
Neumann equation 

dp 



dt 



(34) 



Note that this equation without any explicit time dependence (dp/dt) exp i ic i t leads to an 
unitary time evolution, see also H and Sec. |TV|. 

In this picture the time dependence of the creation and annihilation operators turns out 
to be 



(35) 



However, this trivial time dependence gives rise to the possibility of parametric resonance 
which enhances the chances to verify the effect of particle creation experimentally. Further- 
on we will denote the initial creation and annihilation operators by a M (0) = a M . Note that 
in this picture the particle number operator = d^d^ is time independent for all modes. 



E. Rotating Wave Approximation 



In the interaction picture the time-evolution operator is given by 



U(T,0) = %exp 



-i I [Hf(t)+H]f(t))dt 



(36) 
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where % denotes time-ordering. If the interaction Hamiltonian Hj leads to small correc- 
tions, the usual procedure is to apply perturbation theory via expanding the exponential. 
Evidently, this would imply that U is close to the identity. On the other hand, in order to 
make an experimental verification of quantum radiation feasible, the time-evolution opera- 
tor should deviate significantly from the identity. Therefore a different approximation needs 
to be found. For the case of parametric resonance this may be accomplished via applying 
the rotating wave approximation (RWA), see e.g. |5], [2(], ^TJ. Within this scenario the left 
boundary performs harmonic oscillations obeying a(t) = ciq + e(b — ao) sin(a;i) with the di- 
mensionless amplitude [|38[ | e C 1 and the external vibration frequency to during the time 
interval [0, T\. This also implies an oscillating time-dependence of the frequency deviation 
and coupling matrix 

on . 



»dar 



- a )esin(u;t) + O (e 2 ) , 



M^{t) = m^{b - a )ujecos(ujt) + O (e 2 ) 



(37) 
(38) 



see also Sec. [II Q Together with the trivial time-dependence of the ladder operators (f35|) in 
the interaction picture this enables us to perform the RWA. Expanding the time evolution 
operator (|36D into an infinite series one can treat the time-ordering terms in the following 
way: In analogy to 



% 



+H I (t 1 )H I (t 2 ) 



(39) 



one can rewrite all these terms to yield a multiple product of Hamiltonians without time- 
ordering and terms involving commutators with Heaviside step functions. In the resonance 
case, i.e., when uj = 2Q ^ the terms with commutators can yield nothing but strongly oscil- 
lating integrands ||, as can also be seen by Fourier-expanding the involved Hamiltonians. 
With the duration of the perturbation being sufficiently long, i.e., with ujT ^> 1, the con- 



tribution of these terms to (pq) will be comparably small. As a consequence, time-ordering 
can be neglected to all orders within the RWA. 

The remaining integrals - without time-ordering - factorize and can be resummated to 
yield an effective time evolution operator 



?7 eff (T,0) R ^ A exp {-iH'jT 



(40) 
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where the effective interaction Hamiltonian 

£/ ffT RWA f T \Hf {t ) + HY(t)\dt (41) 
Jo 1 J 

still has to be calculated. The above time integration involves many oscillating terms. Again, 
with the duration of the vibration being sufficiently long ujT ^> 1, i.e., after many oscillations, 
the time-integrated interaction Hamiltonian can be approximated in the following way: Since 
the time average of purely oscillating terms is rather small compared to that of constant 
contributions we may neglect the former ones. As a result, in the series expansion of the 
time-evolution operator U(T, 0) only those terms where the oscillation of the ladder operators 
is compensated by the external vibrations - represented by AQ^(t) and M^ u (t), respectively 
- will be kept. Strictly speaking, in the above equation terms of O (e 7 (u;T) J ) are neglected 
by the RWA if / > J holds. The terms with J = K - i.e., exactly the terms in which the 
oscillations of the creation and annihilation operators (|35|) are compensated by the external 
time dependence Afl^(t) and M^it)] - will be kept. (Note that terms with J < K do not 
occur.) 

The general squeezing interaction Hamiltonian reads 

1 dQ, 

+a^(t)al(t) + al(t)a,(t)] + O (e 2 ) . (42) 

Accordingly, within the RWA only the terms fulfilling the squeezing resonance condition, 
see also e.g. @, 0, \Q |2§ 

u = 2^1 (43) 

will be kept. In general \x can also be a right-dominated mode, but note that in this case 
the effective squeezing Hamiltonian would be of (D(rj 2 ), since according to fl24 ) the right- 
dominated eigenfrequencies do not depend on a(t) up to 0(rj). Therefore we will restrict 
ourselves to left-dominated modes /x and among those in particular to the lowest one, i.e., 
as commonly done we consider the case of fundamental resonance 



2n° lu = 2Q° L . (44) 



From now on this mode will be abbreviated by the index L = (1,1,1,/) throughout this 
publication. Consequently, by virtue of 

i^ T sinM)e^ R ^ A ±^ (45) 
13 



an effective squeezing Hamiltonian can be derived 

H* s = * [(4) 2 - («l) 2 ] , (46) 

where £ is given by 

i /Cpo\ 2 

H^GI) • (47) 

Obviously ffjj is a generator for a squeezing operator for the mode L with £ being the 
squeezing parameter. 

The same procedure can be applied for the velocity interaction Hamiltonian which reads 
in the vibration case 



Hi = \/ "V 6 ^ 6 _ fl o) cos(^t) x 

\M> V " 

-a M (t)at (t) - a M (t)a,(t)] + (e 2 ) . (48) 

However, the occurrence of inter-mode couplings now results in a different resonance condi- 
tion (see also |5|, |22| ) 



K±Q°\. (49) 



Depending on the frequency spectrum of the cavity under consideration this resonance con- 
dition might be fulfilled by several pairs /j,v, but here we will assume for simplicity that only 
one such pair exists. Though in any case via 

i [ T cos(ut)e ±wt dt R ^ A - (50) 
T Jo 2 

an effective velocity Hamiltonian can be derived, two major distinctions should be made. 

a. © coupling uj = + Q®- m this case one yields an effective velocity Hamiltonian 



given by 




H es = i ( \lffi> ~ y fl§ ) rn 12 eu(b - a ) x 

a[a\ — aia 2 ) , (51) 



which is a non-diagonal multi-mode squeezing Hamiltonian. Note that if one wants to 
fulfill squeezing and velocity resonance conditions simultaneously (2VL Q L = uj = + $7°), 
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the number of possible combinations reduces significantly, since then the velocity resonance 
condition can not be fulfilled by two distinct left-dominated modes. For reasons of brevity 
we do not consider this case here. 

b. coupling u = Vt\ — Q®. Here the resulting effective velocity Hamiltonian does not 
resemble a squeezing but a hopping operator 

Iqo' 




(52) 

This coupling is of special interest since if one does not insist on simultaneously fulfilling both 
resonance conditions - parametric resonance might perhaps still be induced by lower external 
frequencies u = — < 2Q° L whose generation would be simpler in an experiment, see 
also Sec. [VQ . In the case of simultaneously fulfilling both conditions several combinations 
may arise. 

• The frequencies both belong to either right- or left-dominated modes. In it 
has been shown that for an ideal cavity with special dimensions, e.g. a cubic one, a 
strong inter-mode coupling can occur. In that case a much smaller exponential particle 
creation rate has been found. Therefore in view of an experimental verification this 
case is counterproductive and not considered here. Instead we propose a cavity with 
transcendental ratios of the dimensions such that there is no resonant inter-mode 
coupling of similar-dominated modes. In any case such a coupling would certainly 
require large quantum numbers of the involved modes. 

• The frequency represents a right-dominated mode and £7° some left-dominated 
mode, respectively. The lowest possible right-dominated frequency Q 2 = would 
then be obtained when = Q° L . As an example, this case will be considered here. 
We want to stress that the used methods are nevertheless applicable to any possible 
combination of couplings. 

Note that the situation would be completely different in 1+1 space-time dimensions where 
- due to the equidistant spectrum - the velocity term always contributes, see e.g. |5|, [24], 



25], [27], |28|, p9] . The coupling right-dominated mode fulfilling fl nx ,n y ,n z ,r = 3fi° wm 



further-on be denoted with the index R = (n x ,n y ,n z ,r). Accordingly, in our considerations 
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the velocity Hamiltonian reads 

H Is = { X (a[a R - a L a R ) (53) 

with 




* = 4 £ "MVi? + vS) (6 - ° )mi '' i (54) 

being the velocity parameter of the system. Since x — ® ( e ^°L m L,R) = O {eVt^rj) [see also 



(29)] it follows that x/£ = (9 (r?) 1 in the limiting case of a nearly perfectly reflecting 



mirror [[39]. This hopping operator is consistent with the visual picture of a semitransparent 
mirror. 



F. Response Theory 



We assume our system to be initially in a state of thermal equilibrium that can be 
described by the canonical ensemble 

exp (-/3H 



p(t = 0) = po 



(55) 



Trjexp {-(3H^ 

with j3 denoting the initial inverse temperature. The expectation value of an explicitly 
time-independent operator Y at time t = T is given by 

(Y(T)) = Tr {fp(T)} 

= Tr J Y% exp 



T 



PoV ex P 



ij I Hj(ti)dti ) x 

T 

/, / H I (t 2 )dt 2 



(56) 



where 7^ denotes the anti-chronological operator (anti-time ordering) - for a more involved 
discussion see e.g. [0. 

If the interaction Hamiltonian represented a small correction, one could expand the time- 
evolution operator into a perturbation series yielding a series expansion for (P(T)). However, 
for the resonance case this procedure is not justified: In the rotating wave approximation in 



Sec. [II E| the interaction Hamiltonian simplified to 



H, 



cff 



it v» 



<^ L aR — did R 



(57) 
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which implies for the time-evolution operator 

U(T, 0) R ^ A exp (-iH 1 ^) . (58) 

The whole expression for computing the expectation value of an operator now becomes much 
simpler 

(Y(T)) R ^ A Tr{ye Xp (-i(4^ + ^)r)p x 

exp (+i(4ff + H^)T) } , (59) 

but since the correction H^T is not small in the case of interest above expression is still 
not practical for applying perturbation theory. Exploiting the smallness of the velocity 
Hamiltonian it will prove useful to separate the two Hamiltonians. This can be achieved 
with the ansatz 

exp (-i(Hi + H^)t) = exp (-iffjr) a(r) (60) 

with a being an auxiliary operator. Differentiation with respect to r yields a differential 
equation that can be solved for a using the initial condition <r(0) = 1. Introducing the 
parameter ordering 7^. in analogy to time ordering {%) the solution for a can be cast into 
the form 



exp I -% I Hes(r)dT 







(61) 



Here the squeezed effective velocity Hamiltonian has been introduced 

7&(r) = exp (+ifljr) flJJ exp (-iffjr) (62) 

which is now dependent on the parameter r. Further-on we shall denote squeezed operators 
by calligraphy letters. By inserting above equations into the expectation value fl5*9"| ) one 
yields 

(Y(T)) R ^ A Tr !^(T)f T exp ^ £ H^dr^j x 

pot? exp (+i J n^(r 2 )dr 2 ^y (63) 
Please note that in this representation also the observables are squeezed 

y{T) = exp (-Htf&r) Y exp (-iff&r) (64) 
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but here using the physical perturbation time T. We will refer to this picture as the squeezing 
interaction picture. Unfortunately the parameter ordering is reintroduced by this procedure 
but as the advantage of these manipulations we are now able to expand the expectation 
value (Y(T)) into a perturbation series with powers of TC^g- Keeping only terms to second 
order one finds 



(Y(T)) = Tr{^(T)p, 
+ Tr \y{T) 



driH^in] 



+ Tr iy(T) / dnH^Wpo / dr 2 H v e5 (r 2 



-Tr^y(T)T T j dnHUn) j dr 2 H^(r 2 )p, 
~Tr h(T)p tJ j dnWSrCn) j dr 2 n^(r 2 ) 



+ O ( (H^) 3 



(65) 



which is now a practical expression for calculating expectation values. 



III. THE QUADRATIC RESPONSE 
A. Squeezing 

According to the results of Sec. [II F| in the squeezing interaction picture both particle 
number operator and the effective velocity Hamiltonian have to be squeezed. The squeezing 
operator 

S{t) = exp f+iH^r) = exp [(d L ) 2 - (a[) 2 ] r) (66) 
implies the following well-known transformation rules (see e.g. ||31|| ) 

b L (r) = S{T)a L S\r) 

= d L cosh(2£r) + d[ sinh(2£r) , (67) 
8£(t) = ^(r)ot^t(r) 

= a[ cosh(2£r) + d L sinh(2£r) , (68) 

which can also be envisaged as a Bogoliubov transformation of the ladder operators. Due 
to the commutation relations (|31|) other modes than the fundamental resonance mode L 
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are not affected by squeezing. Inserting the above expressions into the effective velocity 
Hamiltonian fl53|) one can easily derive its squeezed counterpart 



iXCOsh(2£r) 
+iX sinh(2£r 



(69) 



Note that the squeezed effective velocity Hamiltonian is now dependent on the parameter r 
and that it is still of O (x) = @ iv) which justifies a perturbation-like treatment. 

The same can be done for the particle number operators where again only the fundamental 
resonance mode N L = a\di is affected 



S(T)N l S\T) = S(T)a[S\T)S(T)a L S\T) 
[1 + 2sinh 2 (2£T)] a[a L 
+isinh(4£T) [(a[) 2 + (a L ) 2 



+ sinh 2 (2£T) . 

For brevity we will denote the hyperbolic functions by 

C{T) = cosh(2£T) 
S(T) = sinh(2£T) 



(70) 



(71) 



throughout this paper. 



B. Expectation Values 

Since we are mainly interested in the phenomenon of quantum radiation and thus there- 
fore in the calculation of the cavity particle content after a time T when the disturbance has 
ended. Due to the dynamical disturbance the system leaves the thermodynamic equilibrium, 
see also 0. The quadratic response of the expectation value of the particle number opera- 
tor can be calculated via substituting Y — > N in equation (|B~5|). However, since the initial 
statistical operator p involves arbitrarily high powers in H , it is practical to rewrite the 
expression obtained from (|65|). Utilizing the invariance of the trace under cyclic permutation 
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and the property of time ordering 



A(t),B(t') 



A(t)B(t') + B{t')A(t) 



T 



A(t)B(t') 



A(t)B(t') 



(72) 



the expectation value of interest can be cast into the more convenient form 



(N(T)} = Tr{^(T)p } 



+ Tr 



+ Tr ^ p / dT 2 H eS (r 2 ) 



Po 



Af(T), / dnH^in) 



+ 



% / dnH^in) / dr 2 n^ s (r 2 ),Af(T) 



O ( (Hl s f 



(73) 



This form is now suitable for evaluating the traces since all commutators only concern a finite 
number of creation or annihilation operators. Thus the quadratic response can be brought 
in relation with the initial particle contents of the cavity (Bose-Einstein distributions) 

1 



N° L 



Tr <^ d[d L p 



N r = Trja^a^po} 



e /3Q« _ x 
1 

e^ n °R - 1 



(74) 
(75) 



where (3 stands for the initial inverse temperature of the system. These mean occupation 
numbers incorporate the whole temperature dependence of the quantum radiation - as long 
as the back-reaction of the field on the moving mirror can be neglected. As we shall show 
later, the lowest order term (rj = 0) is in agreement with the results of an ideal cavity, as was 
considered for example in || 0]. Also, since 1% S contains only odd powers of creation and 
annihilation operators for a single mode, the linear response vanishes. Generally, every trace 
involving an odd power of ladder operators vanishes and J\ as well as p do only contain 
even powers. Hence the last trace in Eq. ( [73] ) constitutes the quadratic answer. In contrast 
to an ideal cavity the terms with time ordering are here especially important since they will 
be found to produce leading order terms. 
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C. Particle Creation 



Using the squeezed operators fl6"9|) and (|70|) it is now straightforward to compute the 



(N L (T)) = S 2 (T) 

+ [1 + 2S 2 {T)] N° L 

Y 2 

+ 4tt [3C 2 (T) - 2C(T) - 1 - 2£TS(2T)1 
+ ^ [4C 2 (T) - 2C(T) - 2 - 4£TS(2T)] iV° 
+ ^ [2C 2 (T) - 2C(T)] N° R 

+0 (r] 3 ) . (76) 



As was anticipated, the lowest order term S 2 (T) + [1 + 2S 2 (T)]N® is in agreement with the 
results obtained in || for an ideal cavity. The linear response (in if) vanishes. It might be of 
interest that the leading terms TS(2T) in the quadratic answer stem from the time-ordering 
which is therefore very important. One can see that at long disturbance times T these 
leading terms show the failure of the quadratic approximation since the particle number 
would become negative at some point. This is due to the fact that ( |73|) is a perturbation 
series in J Q T H^(t)dt which will always become large at some time T. This problem can only 
be solved by including all orders in 7], see also Sec. |V|. 

Of course ( ]75| ) can also be applied to the corresponding coupling right-dominated mode 
(whose particle number operator is invariant under squeezing) where one finds 



Again the linear answer is vanishing. For 77 = there would not be any created particles in 
the reservoir due to the dynamical Casimir effect corresponding to a perfect internal mirror. 



commutators and the traces in the expectation value ( f73[ ) 



As a result one finds for particles 



in the fundamental resonance mode L 



(N R (T)) 




(77) 
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It is remarkable that the coefficient of N® in (Nr(T)) equals the coefficient of N R in 
(Nl(T)). As we shall see in Sec. |V B| , this feature is preserved to all orders in r\. 



IV. THE MASTER EQUATION APPROACH 

In this section it is our aim to derive the associated master equation for an effective 
statistical operator accounting for the left leaky sub-cavity or left-dominated modes, respec- 
tively. So far (3+1) dimensional vibrating leaky cavities have only been treated in different 
setups - see e.g. |3(J - where the vibrating mirror is understood as a (quantized) harmonic 
oscillator coupled to the cavity field (the reservoir) or with master equations adequate rather 
for stationary systems - see e.g. fl2|]. It was assumed in [12 that these master equations 



could also be applied when one of the boundaries was moving. The possibility of limitations 
to that procedure as well as the need for a rigorous derivation of the master equation for 
resonantly excited systems have already been expressed in ||12|| . We want to derive such 
an equation starting from first principles. As a test we will also solve the obtained master 
equation and recalculate the quadratic answer for the left mode particles to compare with 
the previous results of Sec. [Ill C[ To obtain a master equation we will closely follow the 



derivation given in [31 



A. Derivation of a Master equation 

Throughout this section we will deploy the squeezing interaction picture where not only 
the time dependence induced by H but also the dependence resulting from H^ s is determin- 
ing the operator time evolution has already been proposed in Sec. |11F[ In this picture the 
time evolution of the statistical operator is governed by a modified von Neumann equation 

^ = -i [*&(*), m] = -&{Mt) • (78) 

The above equation defines the action of the Liouvillian super operator £ on p (see also Ref. 



32|). By defining the projection super operator *P via 

^Y = Pr (0)Tt r \y) (79) 



for all observables Y, where Tr^ means taking the trace solely over the right dominated 
modes we can introduce a reduced density operator accounting for the left-dominated modes 
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only 

p L (t) = TV R {p(t)} . (80) 

Combining above equations it can be shown that the dynamics of the full statistical 
operator p is governed by the Zwanzig master equation 

-z$£(t)it(t,0)(l-$)/3(0) 

-i$£(#p(f) (81) 

where 

t') = exp - $) J £{t")dt"^J (82) 

is the reduced time-evolution super operator. The Zwanzig master equation is exact to all 
orders in r\ but usually too complicated to be solved. However, assuming an initial thermal 
equilibrium state and taking into account that initially our system and reservoir do not 
interact (no correlations) it can be simplified considerably: 

1. In analogy to the argumentation concerning the vanishing of the linear response in 



Sec. [Ill B| it follows that 

Tr R {n^p R (0)} = (83) 

since Tt^s contains only odd and Pr(0) only even powers of the creation and annihila- 
tion operators for the mode R. This can equivalently be written as 

$£(i)$p(t) = . (84) 

2. In our setup the initially stationary system (stationary walls) does not permit inter- 
actions between system and reservoir, since both M^(to) and AQ 2 (to) will vanish. 
Consequently, assuming thermal equilibrium, system and reservoir initially constitute 
independent subsystems which cannot be correlated, i.e., the initial statistical operator 
of the cavity modes factorizes 

p = p(0)=p L (0)®p R (0), (85) 

hence one finds (with Tr R {p R } = 1) 

(l-$)p(0) = 0. (86) 
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These assumptions yield a simplified Zwanzig master equation 



app(*) 

dt 



$£(t) f il(t,t')£(t')fpp(t')dt' (87) 



which is exact but still too complicated to be solved. 

In order to gain a solvable equation we will apply further approximations: 

a. Born approximation 

Since £ = O (rf) one can approximate the reduced time-evolution operator via il(t, i') 
1 + O (rf). This neglects terms of O (rf) if inserted into (|57|) and yields 



$£>) / £(0$P(t')dt' + (r? 3 ) . (88) 
Jo 



By employing the reduced density operator ph(t) = TrR{p(t)} one can equivalently 
write 

° l>L(t) = -Tr R {2(t) f Z(t')p R (V)p L (t')dt 



dt 

+0 (r/ 3 ) . (89) 

This equation governs the time evolution of the effective statistical operator p L ac- 
counting for the left cavity. 

b. Markov approximation 

The retardation in equation (^J), i.e., the occurrence of PL(t f ), complicates the cal- 
culation of p^(t). Iterative application of ( |S9"D implies that Pl(^') = Pl(^) + @{v 2 )- 
Accordingly, we apply the Markov approximation, which is also known as short mem- 
ory approximation, simply by replacing pi(t') — * Pi(t) on the right hand side. Since 
£ = O (rf) we thereby neglect terms of O (r/ 4 ) and yield the Born-Markov master 
equation 

= -Tr R y*2(t)£(t>)p R (0)p L (t)dt^ 

+0 (r) 3 ) , (90) 

thus having maintained the level of accuracy. 
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Using the definition of the Liouville operator £ in Eq. (|78|) one can equivalently write 
dfa(t) 



dt 



+Tr R {H* s (t)p R (0)PL(t) / W 



-Tr„<{W£(f) j( 7^(0^(0)^(0 



+h.c. + (r/ 3 ) . 



(91) 



Finally, having evaluated both traces and after having performed the t' integrations with 
the aid of fl69|), one yields the following master equation 



dt 



fi(t) 2a) L p L {t)a L - a L a\p L {t) - p L (t)a L a[ 
+f 2 {t) 2d L p L (t)d !( L - a ] L a L p L (t) - p L {t)a} L a L 
+fs(t) a/ L p L (t)al + a L p L {t)a L 

-h(t) [(a L fp L {t)+p L (t)(a\f 
+0 (r] 3 ) 



(92) 



where the functions /j(i) are given by 



hit) = ^S(t){C(t)(2N R + l)-N R -l} 



fs(t) 

hit) 



^S(t){C(t)(2N R +l)-N° R } 

2i 
2i 

x 2 



{C 2 (t)+S 2 (t)-C(t)}(2N R + l), 
{[C 2 (t)+S 2 (t)-C(t)]N R + C 2 (t)-C(t)} 



f&) = 7^{[C 2 (t)+S 2 (t)-C(t)]N R + S 2 (t)} 



(93) 



Via averaging over the degrees of freedom of the reservoir and by applying the Born-Markov 
approximation we have now rigorously derived a differential equation for an effective statis- 
tical operator pz,(£) accounting for the leaky cavity. This effective statistical operator obeys 
a non-unitary time-evolution (changing entropy). There are several possibilities to check 
the obtained master equation: As the simplest tests one can verify that the time evolution 



preserves ||40|| the hermiticity and the trace of pi,. 
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A better indication for a correct master equation is the fact that if one takes the limit of 
no squeezing, i.e., in this coupling £ — > 0, the resulting equation corresponds to a harmonic 
oscillator coupled to a thermal bath: With 



limA(t) 
Hm/ a (f) 
lim/ i=3i 4 i5 (t) 



x 2 t(K + 1) 





(94) 
(95) 
(96) 



one arrives at a simplified equation 



dt 



Id 



2a} L p L a L - a L a} L p L - p L ciLcJ L 



+1D- 



K + i 



2a L p L a) L - a} L a L p L - pLa[a L 



+0 (r) 3 ) 



(97) 



Apart from the time dependence of the damping coefficient 7^ = 2\ 2 t above equation is 
exactly the well-known master equation for a harmonic oscillator coupled to a thermal bath, 
see e.g. |33] . 

The time-dependence of 7^ is a remnant of the dynamic master equation describing 
the time-dependent system in the unphysical limit £ — > 0. However, in order to have a 
stronger indication for the correctness of our ansatz we want to solve the master equation 
(|92|) explicitly. 



B. Approximate solution of the master equation 



So far we have neglected terms of O (rj 3 ). The functions fi(t) are already of O (rf) which 
makes it possible to maintain the level of accuracy by applying the additional approximation 
phif) ~ Al(0) on the right hand side of equation (p2[), which could also be envisaged as an 
additional Markov approximation. Accordingly, one is now able to yield a solution for p^ 

h(T) = Pl(0) 

+ Fi(T) 2a[p L (0)a L - a L a[p L (0) - p L {0)a L a[ 

+ F 2 (T) [2a L p L {0)a[ - a[a L p L {0) - p L {0)a[a L 

+ F 3 (T) \a[p L {0)a[ + a L p L {0)a L 
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- F 4 (T) [(a[) 2 p L (0) + p L (0)(a L ) 2 

- F 5 (T) [(a L ) 2 p L (0) + p L (0)(aI) 2 
+ 0(rf) 



(98) 



with Fi(T) = fi(t)dt. Given this effective statistical operator for the leaky cavity one 
is now able to calculate the number of created particles in all left-dominated modes. Note 
that for considering right-dominated modes one would have to derive a statistical operator 
for the reservoir. 



C. Particle Creation 



Since we were working in the squeezing interaction picture - where the observables have 
to be squeezed - the expectation value of the particle number operator reads 

(N L (T)} = Tr L {ri L (T)p L (T)} . (99) 

Other left-dominated modes than the fundamental resonance mode L are trivial to solve: 
Due to the commutation relations ([H]) their ladder operators commute with those of the 
resonance mode L. This implies (due to the invariance of the trace under cyclic permuta- 
tions) that all higher order traces must vanish and one just yields the trivial result of their 
initial occupation numbers. Inserting the approximate reduced density operator obtained 
in Eq. (|98|) as well as J\"l into the above equation, one can see immediately that zeroth and 
first order in t] agree with the previous results but showing this for the second order is a bit 
tedious. After some algebra one finally finds complete agreement which the previous result 
found in Eq. (|7q) of Sec. [Ill Q thus giving a strong indication for the validity of our master 
equation within the RWA approach. 



D. Comparison with other results 



In ]T2| the effects of losses are taken into account by a generalized version of the simple 
master equation ansatz 

dp 



dt 



p,H 



+ — [2dpd^ — a^ap — pa? a] . 



(100) 
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However, as we have observed in the previous calculations, this master equation does not 
adequately describe the leaky cavity under consideration: 

• It is restricted to the case where the initial state of the reservoir is just the vacuum 
state and therefore does not include temperature effects. This has been taken into 
account in HT2 . 



In addition, even the master equation for an harmonic oscillator in a thermal bath KB 



dp 

di 



+ 



p,H 
n + 1 



+ -^1d [2a^ pa — aa) p — paa^ 
7d \2apa) — a/ap — pa^a\ 



101) 



cannot be assumed to describe the system correctly. Even if one identifies the Hamil- 
tonian H in above equation with the effective squeezing Hamiltonian H^ s this master 
equation goes along with serious problems since the Markov approximation is not 
justified anymore. This complication reflects the inherent dynamic character of our 
system. As we have shown in Sec. |1V A| the complete master equation resembles above 



equation only in the limit of no squeezing £ — > - see Eq. (|97 ) - and even then with 
a time-dependent damping constant jjj. 



Instead, the complete master equation ( P2"D displays more similarities to one in a squeezed 
thermal bath where one has to replace the parameters by time-dependent functions. Accord- 
ingly, the dynamical system under consideration is described properly only by an explicitly 
time-dependent master equation. 

Potential limitations to Eq. (|100|) have already been anticipated in . 



V. THE NON-PERTURBATIVE APPROACH 



The previous results in Sec. [Ill C| and Sec. |1V C| have not been able to explain the behavior 
of the system in the limit of a long-lasting disturbance. The leading order term TS(2T) in 
( [75p has a negative sign which would lead to negative particle numbers for large disturbance 
times T. This problem can only be solved by including all orders in rj. In this section we 
present a non-perturbative approach within the RWA which enables a convenient calculation 
of expectation values using computer algebra systems fl3| . As a further advantage we want 



28 



to mention that it can in principle be generalized in a straightforward way to the case of 
more than just two coupling modes as was assumed in Sec. |11 E. 



A. Time Evolution 



Application of the RWA in Sec. |IIE| yielded the effective time-evolution operator U e g = 
exp(— iH^gT) with the effective interaction Hamiltonian (|57D . We want to calculate the 
expectation value of an observable Y 



(Y(T)) 



Tr ■[ Ye- i6 °« T p e +i6 °« T } 



Tr t e +iH ^ T Ye- tH ^ T p 



(102) 



In contrast to the previous sections here the full time dependence is shifted back on the 
operator Y. Since Y can be expressed using creation and annihilation operators and due 
to unitarity of the time-evolution operator one just has to find a solution for the full time 
dependence of the ladder operators which is given by 



a a (T) 



„+iH l „T" -iH'„T 

e eff a a e ctl 



(103) 



The above expression requires special care to evaluate, since H^ s is not a pure squeezing 
generator. In our considerations the effective interaction Hamiltonian is time-independent 
which does not necessarily hold in general. To preserve generality we will therefore introduce 
an auxiliary parameter i9 while keeping the time T fixed. This enables us to write 



(104) 



Obviously we are interested in a a (T) = d a (d = 1). To this end we define a 4 dimensional 
column vector (see also @) 

V 4(^) J 



X(0) 



(105) 



Since H^ s does not depend on one finds 

dx 



d-d 



iT 
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(106) 



where A is a number- valued 4 by 4 matrix acting on x. This form can always be reached 
if the effective Hamiltonian is quadratic: The commutation relations ( |31~D lead to a linear 
combination of creation and annihilation operators that can always be written as a number- 
valued matrix A acting on x. Since A is independent of •& the solution is obtained via 



x(&) = exp (ATti) x(0) 



(107) 



and hence 



x(l) = exp (AT) x(0) = U(T)x(Q) . 



(108) 



Thus the whole problem reduces to a calculation of the time-evolution matrix U_(T) = 
exp(y4T). In the present case the structure of H e s in Eq. (|57[) implies a very simple form of 
.4 



.4 



( 2£ X \ 
2( X 
-x 

V o -x o o y 



(109) 



The four eigenvalues of A are given by 

Ai = 
A 2 = 

A 3 = 
A 4 = 



2 - X 2 



(110) 



For reasons of brevity we shall omit the full listing of the matrix U_(T) = exp (AT) - it 
can easily be calculated using some computer algebra system. Note that the exponential 
matrix U_(T) is positive definite for all T and thus will not exhibit the problems associated 



with the extrapolation of the used approximations of sections |TJ and [TV] beyond their range 
of validity. In order to calculate expectation values one just needs the matrix elements 
of U_{T). This becomes evident considering the time evolution of the new annihilation 
and creation operators x(T) = U_(T)x(0), i.e. Xi(T) = ^2j =1 Uij(T)xj(0). Therefore the 
expectation values of particle number operators of the resonance modes Nl(T) = X2{T)x\(T) 
and Nr(T) = X4(T)£ 3 (T) can be calculated simply by insertion of x(T). After evaluation 
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of the remaining traces containing only the initial creation and annihilation operators Xi(0) 
one finds the full response function to be a combination of matrix elements of U_(T) 

(N L (T)) = (U 12 U 21 + U 14 U 23 ) 

+(U U U 22 + U 12 U 21 )(N° L ) 

+(U 13 U M + U U U 2S )(N° R ) (111) 

and 



(N R (T)) = (U 41 U 32 + U U U 43 ) 

+(U i2 U 31 + U 41 U 32 )(N° L ) 
+(U 33 U U + U 43 U U )(N R ) 



(112) 



B. Particle Creation 

With the full knowledge of U_ we are now in a position to state the full response function 

of the particle number operator. Having inserted the matrix elements of U_ into Eq. ( |1 1 1| ) 

one finds after performing some simplifications 

1 



(Nl(T)) 



4 (e - X 2 ) 



x 



2T e+ 



£ cosh 
+£cosh 2TU- v/? 



2 - X 2 



2 - X 2 



X' 



X 1 



-2 X 2 cosh [2T£] -2(e~X 2 )} 



2(e - X 2 
£ cosh 



x 



2T ( + VF? 



+ecosh 2T u - ve-x 2 ) e - 



2 - X 2 



-x 



cosh [2Tf] cosh (ZZV^ 2 - X 2 ) + 1 } 



K 



2(e - x 2 ) 



x 



cosh [2Tf] cosh ^2Ta/£ 2 - x 2 ) - 1 } 



(113) 
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FIG. 4: Particle creation in the fundamental resonance mode Nl for N® = 1000, = 100, £ = 1 
Hz, and X = 0-5 Hz. An exponential growth is found in all cases. 



This result is valid to all orders in x/£ or Vi respectively. To show consistency with the 
results obtained in Sec. [Ill Q and Sec. |IVC| we expanded the above expression around 
x/£ = up to second order and found complete agreement with Eq. (|76|) ! However, even 
for large values of x/£ = 1/2 the quadratic approximation is a rather good one - provided 
that the duration of the disturbance T is not extremely large - as one can see in Fig. [|. 



Doing the same calculations for the corresponding right-dominated mode one finds as a 
result 



(N R (T)) 



1 



4(e - x 2 ) 



£ cosh 
+£ cosh 



r 



2T[£-y/e 



X 



a + Ve^ 2 



-2 X 2 cosh [zrt] -2(e- x 2 ) } 



Nl 



2(e - x 2 ) 

[ X 2 cosh [2T£] [cosh (2Ty/? - X 2 ) ~ l] } 
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FIG. 5: Particle creation in the right resonance mode Nr for N® = 1000, N R = 100, £ = 1 Hz, 
and x = 0-5 Hz. The lowest order result just corresponds to a constant initial particle number. 
Generally the particle creation in the reservoir is a much smaller effect then in the leaky system, 
see also Fig. 0. 



K 



2(£ 2 - X 2 ) 



£ cosh 
+£ cosh 



2TU + ^e^ 



2T U- 



X 



-X 2 cosh[2T£] cosh(2T v / ^ 



V I + 1 



(114) 

where the remarkable agreement of coefficients of iV^, in (Nl) and of N® in (Nr) as was 
already noticed in Sec. [Ill C] is preserved for all orders in 77. These terms fit the classical 
picture of particle transportation through the leaky membrane where the particle flux is 
proportional to the number of particles on the other side. Again, expanding around x/£ = 
up to second order yields exact agreement with ( fT?|) . Accordingly, also outside the leaky 
cavity particles are produced due to the dynamical Casimir effect, see also Fig. [|. 
Note that at least the quadratic answer is necessary to treat particle creation effects outside 
the leaky cavity. It is still valid that finite-temperature corrections will enhance the pure 
vacuum phenomenon of particle production by several orders of magnitude. (For a direct 
comparison see Fig. [7| in Sec. |VIII| .) 
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FIG. 6: Oscillations of the particle number in the resonance modes Nl and Nr for N® = 1000, 
N R = 100, f = 1 Hz and X = H Hz. 

C. Further Remarks 

We have derived a complete solution for the effective interaction Hamiltonian Q57| ) which 
is valid to all orders in x/£ = (9 (77). As an illustration we consider a case outside our initial 
intentions where x/£ a ^ so assumes large values, e.g. x/£ — !• m this case the arguments 
of the hyperbolic functions in ( |1 13|) and (|114 ) will receive an imaginary part. The arising 



imaginary parts of (N(T)) cancel as they have to because N is a physical observable. Thus 
one finds that if the velocity parameter X exceeds the critical value x ^ £ the particle 
occupation number of the resonance modes versus the vibration time will exhibit oscillations. 
Of course, for the case of a nearly perfectly reflecting mirror inside this scenario is completely 
unrealistic since then x/£, = O iv) wm be relatively small. However, this case is not at all 
academic: If the label R stood for a left-dominated mode - which is the case we excluded 
in our considerations so far and whose equivalent for ideal cavities has been considered in 
- xH ma y ver Y well become large, since ttilr would then be of O (1). 

Similar oscillations of the particle number were also found in the case of strong inter-mode 
coupling in an ideal cavity |7J. 

Note that with x/£ — ® (v) this also leads to an upper bound for the internal mirror 
transmittance f] above which (corresponding to a highly transparent mirror) one finds os- 
cillations that correspond to inter-mode coupling rather than to system-reservoir coupling. 
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From another perspective this phenomenon could also be envisaged as follows: Starting with 
an ideal cavity whose original dimensions do not permit inter-mode coupling one can insert 
a highly transparent mirror (rj ^> 1). This mirror in turn detunes the ideal cavity in such a 
way that it now permits inter-mode couplings as well. 

It is remarkable that in Fig. |6] the phase of the two modes is shifted: When (Nl(T)) 
is at its maximum, then (Nr(T)) is at its minimum and vice versa. This fits nicely with 
the picture of mode hopping mediated by the inter-mode coupling \. One even observes a 
decrease of the particle number in the L-mode for small times. When defining an effective 
temperature [gj this would correspond to an effective cooling of the L-mode. An extreme case 
of this consideration would be the limit of no squeezing, i.e. £ = 0. This would correspond 
to the possibility (see also Sec. [II E| ) of not fulfilling the squeezing but the velocity resonance 
condition. Performing the limit £ — > everywhere in equations (|113|) and ( 114 ) one would 
find pure oscillations of the particle numbers and no exponential growth at all. This case 
is therefore counterproductive for an experimental verification. Note however, that this is 
different for the case of ©-coupling. The consistency with the earlier results leads to the 
conclusion that our approach was justified and the full response function should describe 
the rate of particle production correctly within the RWA. 

Please note that the described procedure also holds for more than just two coupling 
modes: if one has e.g. n modes fulfilling the resonance conditions given in Sec. |11 El the 
formalism still holds and one will have to define a 2n dimensional vector x. Of course then 
n creation and n annihilation operators of these resonance modes will be contained in the 
Hamiltonian and therefore also A as well as U_ will be 2n by 2n matrices. The calculations 
will simply become more involved but can certainly be performed, e.g. by means of computer 
algebra systems. 



VI. DETUNING 



So far we did assume an exact fulfillment of the resonance conditions, i.e. the vibration 
of the left cavity wall did match exactly twice the fundamental resonance frequency Q° L . 
However, in real situations one will of course have to deal with deviations from this desired 
external frequency since it will not be possible to match it with arbitrary precision. In addi- 
tion, the back-reaction of the created quanta might cause the external vibration frequency 
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to change. 

Consequently, we will now discuss the detuned situation - where to assumes slightly off- 
resonant values. For a review of detuning effects see e.g. J7], [8], [12|]. It has been shown in the 
literature that there exist threshold values for the detuning, above which the exponential 
creation of particles disappears. 

Unfortunately the RWA used in our previous considerations cannot simply be generalized 
to this situation. For a slight deviation from the resonance conditions in equations (fHj) and 



( fi9f) the terms with time-ordering - see section [II E| - are no longer negligible in this way. 

We will consider slightly off-resonant situations, where the external vibration frequency 
does not match the fundamental resonance exactly 



CO 



2Q° L (l + 5) 



(115) 



where 5 denotes a small (dimensionless) deviation 5< 1. However, if one considers such a 
variance it is only consequent to include a possible discrepancy of the coupling resonance as 
well, cf. § 



n° R = n° L (3 + A) 



(116) 



where A < 1 denotes the deviation of the coupling right-dominated mode from the 0- 
coupling resonance condition with the fundamental resonance mode. We will adapt the 
multiple scale analysis (MSA) as proposed in [|7| to our scenario of a leaky cavity, see also 
f|. For this purpose we restrict to the results, since the steps in j7|, || can strictly be 
followed - see also appendix |A|. The main difference in these considerations is that we use 
the eigenfunction system of subsection |II C| instead of those of an ideal cavity and that we 
assume the additional deviation (|1 16|) - see also ||. 

In analogy to subsection [V A| one obtains a matrix A governing the time evolution of the 
ladder operators 



x 



a R /V3 
\a R /VsJ 



117) 



The creation of quanta will only be exponential - and thus noticeable - if at least one of the 
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eigenvalues of this matrix 



/ iQ° L 5 



A' 







\ 



V 





2£ -in° L 5 V$x 

-X/VS zfi°(35-A) 

-x/V3 -in^(3«J- A) J 

does have a positive real part. Note that the slight disagreement between the above matrix 
and the one given in Ref. f| is caused by the usage of a different phase (sin instead of cos). 
With the abbreviations 



U = 8f- A X 2 + 12(fl£)*<fA - 2(ft£) 2 A 2 - 20(n° L ) 2 S 2 



and 



V 



+16£ 2 (£ 2 - X 2 ) 

+64(fi° L )Vr + X 2 + (^)V] 
+ (fi° L ) 2 A 2 [8£ 2 + Ax 2 + 52(Q° L ) 2 6 2 + (Qlj 

10 ^ 2 5A[12£ 2 + 8 X 2 + 24{Q° L ) 2 5 2 + 3(ft° ) 2 A 21 



the eigenvalues of the above matrix read (cf. 0, || 

A 



±^\JU±2VV. 



(118) 



119) 



As a consistency check we may set 5 = A = where the eigenvalues reduce to the ones given 
in Eq. ( |110| ). On the other hand, for x — one recovers the usual result of pure squeezing 
in an ideal cavity A + = ±a/4£ 4 — Q 2 L 5 2 . 

Note that in contrast to 0, |8| the inter-mode coupling and thus the parameter x is very 
small x ^ £• This enables us to expand the quantities IA and V into powers of x- The 
condition for a real eigenvalue W + 2\/V > reads 

4£ 2 + fl 2 L SA - 4Q 2 L S 2 



Q 2 L S 2 < 4^-2 



- n 2 L 5 2 + n 2 L (35 - A) ; 



X 



(120) 



Since x is supposed to be small x ^ £ one obtains a significant contribution only if 35 ~ 
A and also in this case merely in the immediate vicinity of the critical value 2£ = Q° L 5. 
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Consequently, the presence of an internal mirror of moderate quality do not drastically 
modify the threshold 

4< H3)' e (l2l) 

for exponential particle creation. However, we would like to emphasize that the shifts of the 
eigenfrequencies of the cavity due to the partly permeable internal wall must be taken into 
account, see also section |IX| below. 



VII. SUMMARY 



We have considered a massless scalar quantum field inside a leaky cavity modeled by 
means of a dispersive mirror. For the case of the lossy cavity vibrating at twice the fun- 
damental resonance frequency we derived an effective Hamiltonian using the rotating wave 
approximation. Within the framework of response theory the magnitude of particle cre- 
ation due to the dynamical Casimir effect was calculated. Furthermore we deduced the 
corresponding master equation via applying the Born-Markov approximation. We found a 
discrepancy to the master equations used so far (see |12j) to describe oscillating leaky cav- 
ities. We also applied a non-perturbative approach for the explicit calculation of the time 
evolution starting from the effective Hamiltonian. All these methods were found to lead to 
consistent results. In addition, the effects of a detuned external vibration frequency need to 
be taken into account. 

It turned out that for the case of moderately low transmission coefficients (or sufficient 
quality factors) the rate of created particles is almost the same as for ideal cavities. The 
squeezing of the fundamental resonance mode as well as the strong enhancement of particle 
production due to the dynamical Casimir effect are preserved in the presence of transparent 
mirrors. 



VIII. DISCUSSION 



In order to illustrate reasonable magnitudes let us specify the relevant parameters: A 
cavity with a typical size of A = 1 cm would have a fundamental resonance frequency of 
Q° L pa 150 GHz i.e., the corresponding coupling right dominated mode must have a frequency 
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of about Q° R = 3Q° L ~ 450 GHz. According to Ref. Jl2| we assume a dimensionless vibration 
amplitude e ~ 10~ 8 . Consequently, in order to create a significant amount of particles one 
would have to sustain the external oscillations over an interval of several milliseconds. At 
room temperature 1/(3 ~ 300 K one finds the initial particle occupation numbers to be 
N® pa 240 and pa 80. Using the above values the squeezing parameter determines to 
f pa 150 Hz. 

As the quality factor Q of a resonator is defined as 

n energy in cavity 

Q = 2ir — , (122) 

energy loss per period 

one finds as a classical estimate yields for our system 

2tt / { 1 Y\ ~ ( 1 



Q= w? =2 *v + wJ ) = °\¥)- (123) 

T denotes the transmission amplitude through the internal dispersive mirror and 7 was 
defined in Sec. II. Assuming a Q-factor of Q ~ 10 8 [O] (and references therein) this would 
imply for the corresponding perturbation parameter 77 = f^/7 = O(10~ 4 ). With these 
values, a reasonable velocity parameter could be given by x ~ 2 mHz. 
The particle content of the leaky cavity is depicted in Fig. 0. 



IX. CONCLUSION 

According to the above considerations it is necessary to vibrate several milliseconds in 
order to produce measurable effects. As already stated, a cavity at finite temperature might 
even be advantageous - provided the cavity is still nearly ideal at its characteristic thermal 
wavelength. However, even after only one millisecond (10 8 periods) a classical estimate 
based on a quality factor of Q = 10 8 would indicate drastic energy losses. On the other 
hand, our calculations based on a complete quantum treatment show that the effects of 
losses are almost negligible compared to the rate of particle creation as long as 77 <C 1. This 



leads to the conclusion that lower cavity quality factors than proposed in |12l , e.g. Q = 10 
[implying 77 = O (10 -3 )] would already completely suffice to justify our approximations [13 . 
Such quality factors are within the reach of the current experimental status. Of course our 
calculations are based on the assumption that the larger cavity - including both the reservoir 
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FIG. 7: Comparison of particle production in the fundamental resonance mode at finite tem- 
perature and without temperature effects. At room temperature (300 K) the initial occupation 
numbers result as N® = 240 and = 80. Accordingly, squeezing and velocity parameters are 
given by £ = 150 Hz and x = 1 mHz. At room temperature the particle number reaches significant 
values much faster. 




FIG. 8: Sketch of a vibrating cavity enclosed by a larger one. This configuration may facilitate 
the experimental verification of the dynamical Casimir effect inside the smaller cavity. 

and the leaky cavity - is perfectly conducting. The error made by this presumption is of 
O (Q~ 2 ) and therefore certainly negligible. Consequently, the experimental verification of the 
dynamical Casimir effect could be facilitated by a configuration where the vibrating cavity is 



enclosed by a slightly larger one as is demonstrated in Fig. H. A further important result JT3 
of our investigations is the shift of the cavity eigenfrequencies ( p4"[ ) of O (rj) = O (10~ 3 ) which 
needs to be taken into account in order to make an experimental observation of quantum 
radiation feasible. 
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X. OUTLOOK 



A. Multi-Mode Coupling 

In section [11 E| we assumed that only one of the right-dominated modes fulfills the res- 



onance condition for the velocity Hamiltonian, i.e., exactly two modes are coupled. If, 
however, the reservoir cavity becomes larger, the spacing between different levels of its spec- 
trum decreases so that eventually more than just one right-dominated mode begin to couple 

- at least within the range of detuning. In this case the effective velocity Hamiltonian 
would constitute a sum of single two-mode coupling Hamiltonians as the one in Eq. fl53p 

- but accounting for different right modes Ri, R2, etc. As we have observed in Fig. ^ in 



Sec. [VB| , the quadratic answer is completely sufficient for reasonable values of x/C Insert- 
ing the aforementioned sum of Hamiltonians into the quadratic answer one observes that 
the mixing terms vanish. Since the effective velocity Hamiltonian only contains odd powers 
of the creation and annihilation operators per mode, one can only obtain a non-vanishing 
trace if it involves two operators of the same (right) mode. Therefore the quadratic answer 
also decomposes into a sum of contributions each accounting for one right-dominated mode. 
Hence we expect the general structure of (Nl(T)) to persist - just substitute \ 2 ~^ J2iXi i n 
the leading contributions [[41]. In order to ensure the applicability of the perturbative treat- 
ment the number n of coupling right-dominated modes has to be small enough to satisfy 

Typically the spacing between two neighboring right-dominated modes is of O (1/L), i.e., 
the inverse of the characteristic length of the reservoir. On the other hand, the width of 
the resonance peak is of O (1/T). Consequently, if the duration of the disturbance T (1 ms) 
exceeds the characteristic length of the reservoir (which is the case for L < 10 5 m) then n is 
certainly small enough. 



B. Electromagnetic Field 

So far, we have considered a noninteracting, massless, and neutral scalar field. The next 
step could be to extend the calculations to the electromagnetic field. In this case several 
new difficulties arise: 
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1. The boundary conditions cannot just simply be described by Dirichlet (or Neumann) 
conditions. Especially for moving walls their form will be more complicated due to 
Ampere's law (mixing of E and B). 

2. As the electromagnetic field is a gauge theory, one has to eliminate the unphysical 
degrees of freedom in order to quantize it. Again, for dynamic external conditions this 
requires special care, see e.g. p2"f . 

3. The different polarizations of photons need to be taken into account which are of 
special interest concerning the fulfillment of the resonance conditions. 



According to [35] the eigenmodes of the stationary cavity can be divided into TE and TM 
modes. For several cavities (rectangular, cylindrical, spherical) the eigenfrequencies are well- 
known. This enables one to determine the squeezing part of the interaction Hamiltonian. 

In order to deduce the velocity Hamiltonian it will be necessary to find an appropriate 
model for the dispersive mirror. This can be achieved by using a thin dielectric slab with a 
high permittivity: e(x) = 1 +7<5(x). As has been shown for a stationary system in [FJJ] this 
leads to a similar eigenvalue equation as (|23|). 



For the detection of the created field quanta some detecting device will have to be placed 
inside the cavity, e.g. an atom. However, the detector will always influence the created field 
as well. A simple approach for the modeling of a two-level system has been provided in 



28fl . In addition, the non-adiabatic parametric modulation of the atomic Lamb shift 



as has been considered in |37[] - must be taken into account, since it will cause excitations 
of the atom as well. 

Note that the induced quantum field will also excite the internal degrees of freedom of 
the cavity mirrors - an alternate description of losses should therefore also take the energy 



dissipation of the losses within the mirrors into account, see e.g. |p6 |. 

Future work combining all these effects is of immense importance regarding experiments 
on quantum radiation using the dynamical Casimir effect. 
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APPENDIX A: MULTIPLE SCALE ANALYSIS 

Starting with the Lagrangian (|3]) it is straightforward to show that the field operator 
fulfills a modified wave equation 

{n + 2V(r;t)}tp(r,t) = 0. (Al) 

If one now follows by introducing ladder operators via the expansion 

V(r,t) = j2^u n (r,t) +h.c, (A2) 

n 

where 

u n (r,t<0) = -*f n (r)e-** t (A3) 

u n (r,t>0) = J2Qk\t)fk(r;t), (A4) 
k 

one can derive a time evolution equation for the coefficients Q1{t). Using the properties ([5]) 
of the eigenfunctions one obtains 

k 

k 

+\ 2 {t)Y,gik9imQt\t) 

k,i 

h 

(A5) 

where X(t) = L x (t) / L x (t) and L x (t) = c — a(t) in our scenario. The antisymmetric coupling 
g m k is defined via 

9mk = L x [ ^(r)f k (r)d 3 r (A6) 

./cavity OL/ x 
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and is therefore related to the geometry factor m m k via g m k = —L x m m k- Note that compared 
to |7| the last line in ( |A5| ) constitutes an additional term, since in our scenario the coupling 
between different modes may depend on the cavity parameters, see also equation (p9|). 
However, this difference is of minor relevance, since all these terms are accompanied by a 
factor of A 2 (t). If one assumes periodic oscillations of the cavity L x (t) = L x [l + esin(fit) + 
ef(t)}, these terms can be neglected if the amplitude e is small. (The auxiliary function f(t) 
is chosen to meet the continuity conditions on L x (t), see also @.) Consequently, one can 
expand (|A5|) in powers of e 1 to yield 

0l n) (*) + (n2) 3 Qi B) (*) = -2nl^L x esm(Qt)Qt\t) 

j 

+2eVt cos(fit) 9kjQt ] (*) 
j 

+eO(f) + 0(e 2 ) . (A7) 

This equation completely resembles the one found in [0. Note however, that we have to 
use the shifted eigenfrequencies and the eigenfunctions for leaky cavities. An approximate 
solution - for a more detailed discussion see @ - can be obtained via introducing a new 
time scale r = et and inserting the formal expansion 

Qt\t) = Qt m {t,r) + eQt ){1 \t,r) + O (e 2 ) (A8) 
with the unknown functions Q^ 1 ^ ^ 1 ^ into equation ( |A7| ). Finally, one has to sort in powers 



of e. To lowest order one finds a free harmonic oscillator which can be solved by 

QjW<°> (t, t) = 4 n) (r)e^ + flW (r)e-^* . (A9) 

The next order terms (proportional to e) yield a driven harmonic oscillator equation for 
Qk w ith the eigenfrequency Q° 

dlQP w + (^) 2 Ql" )(1) = -29 2 r gi" )(0) 



-2n° k L x ^sm(nt)Q^ m 
oL x 



+2fi cos(Qt) Y 9kjd t Qj 



(n)(0) 



(A10) 
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In order to keep the expansion ( |A8|) convergent, above oscillator must not be at resonance. 
Consequently, all terms proportional to exp(±zf2°t) - with k being the particular mode of 
interest - on the right hand side have to cancel. By imposing these conditions for the mode 
k = L and for the coupling mode k = R and inserting the frequency deviations 



n 



2tt° T + h 



I ZT — QOO 



2ft? + ea 



m u L + H = 3tt" L + e/3 . 



(All) 
(A12) 



one finds four linear and coupled evolution equations for the coefficients A^/(t), B^\t) 

couplings and th 
Having applied the 



A r (t), and B^\t). These equations are - apart from the different couplings and the 
additional deviation H - virtually identical with those presented in [ 
modified phase transformations 



4V) 



^iar/2 (n) 



r) 



r/2 b (n) { 



e iar '-o\ 'IT) 



4 n) (r) 

Bt\r) 



e+ 3iar/2 e - if 5r a (n)^ 
e -3iar/2 e+ i(3r b (n)^ 



(A13) 
(A14) 
(A15) 
(A16) 



it is straightforward to rewrite these equations in matrix form. 



d_ 

dr 



6<"> 

(n) 



M 



6<»> 



where 



M 



% 

Ti 



a 71 
i 

+ 2 a 



\ 



-73 






-7s 



72 





--(3a -2(3) 



72 





--(3a -2/3) J 



The quantities 71,72,73 have been introduced for convenience 



7i 



L x dn° L r 1 (nf) 2 

2 dL x ~ 2 



(A17) 



(A18) 
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(2fi° +h)(2n° L + H-h/2) 

72 " 2& L 9LR 

~Wl9lr, (A19) 
(2ft° + h)(2Q° L + h/2) 2 n0 

73 - 2(3^ + H) 9lr « gfl^Lii 

« | . (A20) 
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